Thermal noise of whispering gallery resonators 
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By direct application of the fluctuation-dissipation theorem, we numerically calculate the fun- 
damental dimensional fluctuations of crystalline CaF2 whispering gallery resonators in the case of 
structural damping, and the limit that this noise imposes on the frequency stability of such res- 
onators at both room and cryogenic temperatures. We analyze elasto-optic noise - the effect of 
Brownian dimensional fluctuation on frequency via the strain-dependence of the refractive index - 
a noise term that has so far not been considered for whispering-gallery resonators. We find that 
dimensional fluctuation sets a lower limit of 10"^'' to the Allan deviation for a 10-millimeter-radius 
sphere at 5 K, predominantly via induced fluctuation of the refractive index. 

PACS numbers: 42.60.Da, 05.40.Ja, 06.30.Ft 



I. INTRODUCTION 

The properties of a body in equilibrium at a finite 
temperature undergo continuous fluctuation about their 
mean values. The fluctuation-dissipation theorem of 
Callen and Welton [l| established the relationship be- 
tween these fluctuations and the dissipative properties of 
the system. This fundamental relationship is of great 
practical importance in precision measurements, such 
as the measurement and stabilization of the frequency 
of laser light using high-quality-factor optical cavities. 
Thermal noise sets a limit on the frequency stability of 
optical cavities, and hence on their performance in such 
applications. 

The thermal noise of conventional Gaussian-mode op- 
tical cavities has been analyzed in detail [^Tal, provid- 
ing guidance in the design of such cavities. A relatively- 
new type of high-finesse optical cavity is the whispering- 
gallery resonator ( WGR) |6| , in which light is confined in- 
side an axially-symmetric dielectric object by continuous- 
total-internal reflection. Solid dielectric optical WGRs 
were quickly recognized as potentially useful [7J due to 
their high quality factors and small mode volumes, result- 
ing in an intense confined coherent optical field. The low- 
est cavity loss demonstrated so far resulted in a linewidth 
of about 1 kHz in pure crystalline CaF2 Q , making these 
resonators very attractive candidates for precision laser 
frequency stabilization. Understanding thermal noise in 
such resonators is critical for this application. For a 
WGR, both the fluctuation of the radius and of the re- 
fractive index must be considered. We present a calcula- 
tion, by direct application of the fluctuation-dissipation 
theorem, of the frequency instability due to fundamental 
thermal fluctuation of the dimensions of crystalline CaF2 
resonators. 

Initial estimates of thermal noise in WGRs have been 
made [9|, [lOj , in particular for thermorefractive noise, the 
fluctuation of refractive index due to fundamental tem- 



perature fluctuations. The thermorefractive noise has re- 
cently been measured for crystalline resonators at room 
temperature [ll| . The thermal noise due to fundamental 
dimensional fluctuations at a fixed temperature has been 
less investigated. For laser frequency stabilization, the 
fluctuations on timescales of 0.01 s to hundreds of sec- 
onds are of interest, for applications such as clock lasers 
for optical clocks [iJI or tests of fundamental symmetries 
such as Lorentz invariance [ij| . We flnd that at low tem- 
peratures the thermal noise due to dimensional fluctua- 
tion dominates over the thermal noise due to fundamen- 
tal temperature fluctuation. We find that the dominant 
component of the dimensional-fluctuation thermal noise 
is due to induced changes of the refractive index, via the 
elasto-optical effect. 

To calculate the thermal noise we follow the direct ap- 
proach introduced by Levin 3] for standard two-mirror 
resonators, adapting it to WGRs. We perform the calcu- 
lations numerically, by finite-element analysis, for sphere- 
and disk-shaped resonators. 



II. CALCULATION 

A. Temperature fluctuation 

Our calculations reported in this paper are focused on 
the fundamental dimensional fluctuations that occur at a 
fixed temperature. However, to determine the expected 
dominant thermal noise terms in different regimes and 
to make contact with existing thermal noise results for 
WGRs, we compare our calculated noise levels to thermal 
noise arising from fundamental temperature fluctuations 
at a given mean temperature. These cause WGM fre- 
quency fluctuation via the dependence of refractive index 
on temperature (thermorefractive noise, TR) and ther- 
mal expansion of the resonator (thermoelastic noise) . We 
evaluate the thermorefractive noise of a sphere using the 



expression of Gorodetsky and Grudinin 
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FIG. 1. (a) Simulated mode profile of a CaF2 disk, as sliown 
by tfie electric field intensity (arbitrary units). The radius 
(R), surface radius of curvature (S) and thickness are all 1 
mm. (b) Simulated strain energy density distribution in the 
same disk, (c) Geometries used for disk simulation. 
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where F evaluates to 0.847 for a l-millimetre-radius CaF2 
sphere at low frequency and for the entire temperature 
range of interest, and 7 is the thermal conductivity. 



B. Brownian boundary noise 

The fluctuation-dissipation theorem says that the 
power spectral density of fluctuations of coordinate x of 
a body, at frequency / is 
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where Y{f) = 2-Kifx{f)/F{f) is the susceptibility for 
displacement x{f) under action of the force conjugate to 
it, F{f). Levin Q noted that 
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where Wdiss is the dissipation when an oscillating force of 
magnitude F conjugate to the desired displacement and 



of frequency / is applied. The cycle-averaged dissipation 
can be written as 
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where U is the strain energy in the body at the tem- 
poral maximum of the applied force and <^(/, T) is the 
loss angle of the material. We assume that </> is constant 
with respect to frequency (structural damping) . The con- 
stant loss-angle model is supported by experimental mea- 
surements on ultra-low-expansion material (ULE) Fabry- 
Perot cavities in the frequency range 1 Hz down to 3 niHz 
[3j. We thus have 
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where in ([7]) we have specialized to the case x = R and 
divided by the radius squared to get the spectral density 
of relative radius fluctuations, and the /-dependence of 
U has been dropped (static approximation). Since the 
relative change of the WGR eigenmode frequency v with 
mode radius R is dv/v = dR/R [la], the Allan deviation 
of the relative resonance frequency fluctuations for this 
1// spectrum is independent of integration time, 
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We refer to this type of thermal noise as "Brownian 
boundary" (BB) noise. 

Thus, to calculate the fluctuation of the mode path 
length, we apply an oscillating force conjugate to a 
stretching of this path. This is a force along the circum- 
ferential optical path in the interior of the resonator, with 
a transverse spatial distribution identical to the mode in- 
tensity profile. A circumferential (ring) load is elastically 
equivalent to a radial one JJ:] . As the mode is very close 
to the surface of the resonator, the load can be closely 
approximated by a surface radial stress distribution with 
the profile of the 1-dimensional (radial) integral of the 
mode profile. This is the load that we apply. 

The calculation is carried out in two steps. First, the 
electromagnetic field of a fundamental whispering gallery 
mode (chosen to be at ^ 1.9152 x 10^^ Hz (1565 nm) for 
concreteness) is simulated using the PDE module of the 
finite-element package COMSOL [17], using weak-form 
partial differential equations following the procedure of 
Oxborrow [18,, J^] . The field intensity along the vertical 
midline of a 2D cross-section of the mode is then fit by 
a gaussian. The field intensity distribution in the radial 
direction is also fitted, for use in the elasto-optic calcula- 
tion described in the next subsection. These fitted pro- 
files are the desired information from the electromagnetic 
simulation. 

The second step is to calculate the strain energy stored 
in the resonator when a load with the same profile as 



the electromagnetic mode intensity is applied to the res- 
onator surface. We do this for a resonator with its axis 
along the [111] crystal direction, and in the static approx- 
imation (i.e. material acceleration is neglected), accurate 
for low-frequency noise. The calculation is carried out 
using COMSOL's 3D structural mechanics module. The 
strain energy is integrated over the body to give U and 
the applied load is integrated over the surface to give F 
[20| . Eq. ([8]) then furnishes the Allan deviation. In both 
the electromagnetic and strain-energy simulations, care 
is taken to mesh sufficiently densely in the mode region. 
In the strain energy simulation, deformation patterns ob- 
served corresponded well to the applied load profiles, in- 
dicating that the force magnitudes used were sufficient. 
Numerical values of CaF2 properties used in the calcula- 
tion are listed in Appendix A. 



C. Elasto-optic noise 

The dimensional fluctuation of the resonator also af- 
fects the mode frequency indirectly via the the elasto- 
optic effect, the dependence of the refractive index on 
strain. Relative fluctuations of different points in the 
resonator material correspond to fluctuating strains, and 
the refractive index of the material varies with strain as 
described by its elasto-optic coefficients f2l| . We refer to 
this type of thermal noise as "elasto-optic" (EO) noise. 
This noise term has also recently been considered for 
multi-layer mirror coatings in Fabry- Perot cavities [22l |. 
The effect of the fluctuation of strain component e on the 
Allan deviation can be written as 
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where p is the elasto-optic coefficient {p = gj^ , where 
e is the dielectric constant), and we have used the ap- 
proximation ^ ~ — ^, valid for azimuthal mode index 
> 1. 

As for the BB noise, we calculate a^ using the direct 
application of the FDT. To calculate the fluctuation of 
a given strain component, we apply a load to the res- 
onator conjugate to that strain that is spatially weighted 
according to the intensity of the mode. We evaluate the 
strain energy in the resonator resulting from the applica- 
tion of this load, and then obtain the strain fluctuation 
by using for x in Eq. ^ the dimension corresponding to 
the strain in question. 

In general there are six strain components 
(en, £22, £33, ei2, £13, £23) to be considered. We re- 
strict our consideration to volume fluctuations, 
and we therefore consider only the dilational strain 
dV/V = £11 + £22 + £33. We expect from the elasto-optic 
tensor that the effects of shear strain fluctuations will 
be of similar magnitude. 

We perform the dilational calculation by writing the 
volume of the mode as ~ 27r^r^i? where r is the mode 



toroid effective minor diameter. Thus the relative fluctu- 
ation of the mode volume is 



'^dV/V — <^dR/R + '^<^dr/r 
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o'dB./B. is given by the calculation of the previous section 
{o'dB./R = o'bb)- We calculate adr/r by applying a stress 
distribution S radially directed to the mode center. 
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where the mode center in cylindrical coordinates (p, 9, z) 
is the circle (po,0, 0). Then adr/r is given by 
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For r in the above relation we use the (averaged) gaussian 
1/e intensity radius r — ^WpWz- The dilational elasto- 
optic coefficient is given in terms of the cartesian tensor 
components by (pn -|- 2pi2)/3. Thus from Eq. ([TU)) and 
Eq. (m we have i23 , ] 
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In our calculations we find that <Jdr/r/<^dR/R ranges from 
12 to 160 for the resonator sizes considered, and we there- 
fore neglect the ffist term in Eq. (fT4| . The contribution 
of the (TdR/R. term increases as the WGR size decreases. 



III. SIMPLE ESTIMATES 

Brownian boundary noise - We can make a simple esti- 
mate of the BB noise of a spherical resonator by con- 
sidering the strain energy that results from a uniform 
pressure on its surface. We expect that this will yield a 
noise smaller than the true value for a spherical WGR, as 
in this estimate the radius fluctuations will be averaged 
over the entire spherical surface surface rather than over 
only the mode region. For a sphere subject to a constant 
surface pressure P the strain energy U is 



U = P/W = P 



VP An P^R^ 
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where k = (Cn -I- 2C12) /3 is the bulk modulus for a 
cubic crystal with Cn and C'12 elastic constants. The 
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applied force is F ^ AttR'^P, thus we estimate 



^bb{t) 



For a l-millinietre-radius sphere at 5.5 K, with 
(j}{bK) = 2 X 10-8 m, this gives asB = 3 x IQ-^^ for 
CaF2 (k = 90 GPa), a factor of 2.5 smaller than the 
value determined by more detailed calculations below. 

Elasto-optic noise - We can model the WGM volume as 
a tube of length 27ri? and radius r. We estimate r from 
Vm = 2Tr^ Rr^, where the mode volume Vm is given by 
Kn ~ 3.47r3/2(A/n)^/6i?ii/6^ where A is the optical 
wavelength, 



r ~ 0.335(A/n)7/i2i?5/i2 



We calculate the strain energy that results from a uni- 
form pressure P applied to the cylindrical surface of the 
tube, approximating the WGR material as isotropic. Ap- 
plying the plane-strain condition, the non-zero stresses in 



the tube are ar 



-P. cr. 
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Poisson's ratio, and the resulting strain energy is 

3 



U = n^r^RP^ 
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of the cylinder's fundamental drum mode as representa- 
tive of the bulk material loss and therefore as the limiting 
value achievable for losses. The Q of this mode ranged 
from 10'^ to 5 X 10® between 5 K and 300 K. Using this 
loss angle = 1/Q = (2 x 10-^ 5 x 10"®) (5 K, 300 K 
values), we obtain the Allan deviation values shown in 
Table IH Uncertainty in cj) is the principal uncertainty in 
our calculation. For different values of (j> the BB and EO 
noise terms scale as 0^/^. 

Figure [2] shows the different thermal noise terms for 
a millimetre-sized resonator. At room temperature and 
short time scales, the dominant thermal noise contri- 
bution is the thermorefractive noise, in agreement with 
experiments [l3, \lM- The BB and EO noises vary 
with temperature as [T(j){T)]^^^, while the temperature- 
dependence of the thermorefractive noise is the much 
stronger function T{l/n){dn/dT)y/j{T). Thus the BB 
and EO noise terms become relatively more important 
with decreasing temperature. The EO noise is an order of 
magnitude larger than the BB noise for millimetre-sized 
WGRs, and becomes the dominant noise component at 
low temperatures and long time scales. 

Figure[3]shows the scaling of the calculated noise terms 
with disk radius and surface radius of curvature. We see 
that within the range of values considered, the noise is es- 
sentially independent of surface radius of curvature, and 
the ratio of EO to BB noise scales approximately as i?°'^ 
(19) (fEo oc ~iZ-°®, aBB oc ~i?-^'^). For a sphere the seal 



(18) 



where G is the shear modulus. The applied force is i^ = 
(27r)2ri?P, and we get 
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For a 1-mm-radius resonator, at A = 1.56 /xm and with 
G = 45 GPa, this numerically evaluates at 5.5 K to 
7 x 10^^^, close to the more accurately calculated result 
below. 



IV. RESULTS 

We perform the calculation for sphere and disk radii 
ranging from 0.1 mm to 10 mm, and for disk surface radii 
of curvature ranging from 0.1 mm to 10 mm. The disk 
shapes used are shown in Fig. |TT1 To numerically evalu- 
ate the Allan deviation with equations ([8|) and (fT4|) we 
must know the material loss angle (j). Nawrodt et. al. 
[24| measured the Q of mechanical modes of a pure crys- 
talline CaF2 cylinder as a function of temperature. In 
the cylinder's fundamental drum mode, material motion 
is concentrated near the cylinder axis and away from the 
experimental support. We therefore take the reported Q 



ing is similar {aEO oc ~ i? 



-0.9 
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*). This 



scaling agrees with the simple estimates in the previous 
section, for which aso oc R~'^'^, cbb oc R^^^^. The 
thermorefractive noise depends on mode volume in the 
same way as does the EO noise, and should thus have a 
similar dependence on R. 

Thus we see that lower noise levels are achieved by 
going to lower temperatures and larger resonators, with 



TABLE I. Calculated Brownian boundary and elasto-optic 
noise for CaF2 spheres and disks at 300 K and 5 K. R: radius; 
S: surface radius of curvature. 



R, S 


cbb 




O'EO 




cr_BS 




O'EO 




(mm) 


RT 




RT 




5K 




5K 




Spheres 


















0.1 


2x10' 


-14 


9x10" 


-14 


2x10" 


-15 


8x10" 


-15 


1.0 


8x10- 


-16 


1x10" 


-14 


7x10" 


-17 


1x10" 


-15 


10 


3x10' 


-17 
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FIG. 2. (Color online) Thermal noise for a 1-millimetre-radius 
CaF2 sphere from 5.5 K to 300 K. TR: thermorefractive noise; 
BB: Brownian boundary noise; EO: elasto-optic noise. There 
is a zero in the first-order thermorefractive noise at ~33 K 
due to the zero-crossing of dn/dT. 
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FIG. 3. (Color online) 5.5 K Allan deviation due to Brownian 
boundary noise and elasto-optic noise of a CaF2 disk WGM, 
as a function of disk radius R and surface vertical radius of 
curvature S. BB, R — 1.0 mm: Brownian boundary noise as 
a function of S, for R = 1.0 mm. ih = 2 x 10~*. 



extrapolated noise levels at the 10^^^ level for a lO-mm- 
radius sphere at 25 mK and for a 10-cm-radius sphere at 
1.7 K, where we have used the 5 K value of 6. 



V. CONCLUSION 

We have calculated the thermal noise limit to the fre- 
quency stability of crystalline CaF2 whispering-gallery 
resonators set by dimensional fluctuations, by direct ap- 
plication of the fluctuation-dissipation theorem and us- 
ing the widely-applied structural damping model. We 
have identified a new source of thermal noise in whis- 
pering gallery resonators, due to the elasto-optic effect. 
This noise is smaller than the thermorefractive noise at 
room temperature, so it has not been observed in pre- 
vious room-temperature experiments. At temperatures 
sufficiently below room temperature (e.g. 50 K), the 
elasto-optic noise dominates. The total thermal noise 
is reduced by increasing the resonator size and reducing 
the temperature. WGRs have already been operated at 
cryogenic temperatures [23, l2g| : and therefore an experi- 
mental study of the noise at low temperatures should be 
feasible. One important result of such studies will be a 
determination of the loss angle 4>, which can at present 
only be estimated. For achieving extremely low total 
noise, it would be favourable to operate large resonators 
{R > 50 mm) at temperatures < 100 mK, where a noise 
level < 4 X 10~^^ is estimated from the present work. 



Appendix A: Material properties 



TABLE II. CaF2 material properties. Property values at in- 
termediate temperatures may be found in the cited references. 



Property 




Value 


Ref. 


Refractive index (n) 




1.43 


mi 


Elastic constants: 


Cn 


164 GPa 


[281 




C12 


53.0 GPa 






C44 


33.7 GPa 




Thermal conductivity (7): 


300 K 


9.71 Wm-^K-i 


[2J 




5.5 K 


1200 Wm-iR-^ 




Thermorefractive index 


295 K 


8.4 X lO-'^ 


[271 


{\{l/n){dn/dT)\): 


5.5 K 


1.5 X 10"^ 




Loss angle {4>): 


300 K 


5 X 10"** 


[241 




5.5 K 


2 X 10"* 




Elasto-optic constants: 


pii 


0.039 


[301 




P12 


0.223 






P44 


0.051 




Thermal expansion 


295 K 


1.89 X 10"^ K"^ 


^1 


coefficient (a): 


5.5 K 


6.4 X 10"" K-' 


a 



Extrapolation from 20 K data using T^ dependence. 



Appendix B: Calculation results 
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a. Electromagnetic modeling 



TABLE III. Modeled fundamental electromagnetic modes: 
sphere 



Sphere 


Mode 


Mode 


Mode polar 


radius 


azimuthal 


frequency 


intensity 1/e^ 


(mm) 


index 


( 10" Hz) 


half-width (/im) 


0.1 


559 


1.9164 


4.25 


1 


5706 


1.9152 


13.5 


10 




1.9149 


42.0 



TABLE VII. EO noise calculation for CaF2 spheres. Wp-. 
mode radial 1/e^ intensity half- width, pg: radial position of 
mode center. 



R 


Wz, 


Wp 


Po 




U 




F 


f^dr/r 




(mm) 


(pm) 


(mm) 




(J) 




(N) 


at 5.5 K 




0.1 


4.25, 


1.1 


0.0986 


7.1 


xlO" 


-14 


0.0088 


2.3x10- 


14 


1 


13.5, 


2.5 


0.996 


4.6 


xlO" 


-ii 


0.65 


2.9x10" 


15 


10 


42.0, 


5.0 


9.99 


2.3 


xlO" 


-8 


41 


4.2x10^" 


16 



TABLE IV. Modeled fundamental electromagnetic modes: 
disk 



Disk Vertical radius Mode Mode 


Mode polar 


radius of curvature azimuthal frequency 


intensity 1/e^ 


(mm) (mm) index ( lO" Hz) 


half-width (/im) 



0.1 


0.15 


559 


1.9156 


4.6 


1 


0.15 


5706 


1.9152 


8.2 


10 


0.15 


57320 


1.9151 


14.1 


1 


0.10 


5707 


1.9151 


6.7 


1 


1 


5707 


1.9151 


11.3 


1 


10 


5708 


1.9155 


13.0 



b. Strain energy calculation 



TABLE V. BB noise calculation for CaF2 spheres. R: sphere 
radius, Wz'. mode 1/e intensity half- width, U: strain energy 
10® N/m^ peak load, F: total force, (JciR/R: Allan deviation 
for = 2 X 10"*. 



R 

(mm) 



0.1 
1 
10 



Wz 



4.25 
13.5 
42.0 



U 
(J) 



3.2x10"" 
4.4x10"" 
5.3x10"® 



F 

(N) 



0.0047 
0.15 

4.7 



(^dR/R 

at 5.5 K 
2.0x10"''' 
7.2x10"'^ 
2.6x10"^*^ 



TABLE VI. BB noise calculation for CaF2 disks. R: disk 
radius, S: surface radius of curvature. 



TABLE VIII. EO noise calculation for CaF2 disks. 



R, S Wz, Wp po U F 

(mm) (/im) (mm) (J) (N) 



^ dv I r 

at 5.5 K 



0.1, 0.15 


4.6, 


1.1 


0.986 


8.8x10"'"* 


0.0094 


2.3x10" 


-14 


1, 0.15 


8.2, 


2.5 


0.997 


1.2x10"" 


0.39 


3.3x10" 


-15 


10, 0.15 


14.1, 


5.3 


0.999 


1.5x10"® 


14 


5.2x10- 


-16 


1, 0.1 


6.7, 


2.5 


0.997 


6.7x10"'^ 


0.31 


3.3x10" 


-15 


1, 1 


11.3, 


2.5 


0.997 


2.9x10"" 


0.53 


3.1x10" 


-15 


1, 10 


13.0, 


2.5 


0.997 


4.2x10"" 


0.61 


3.0x10" 


-15 



-ft, o 

(mm) 



Wz 

(^m) 



U 

(J) 



F 

(N) 



O'dR/R 

at 5.5 K 



0.1, 0.15 


4.6 


3.8x10" 


-la 


0.0051 


2.0x10" 


-15 


1, 0.15 


8.2 


1.9x10" 


-11 


0.091 


7.7x10" 


-17 


10, 0.15 


14.1 


1.8x10" 


-9 


1.6 


4.5x10" 


-18 


1, 0.1 


6.7 


1.3x10" 


-11 


0.075 


7.9x10" 


-17 


1, 1 


11.3 


3.3x10" 


-11 


0.13 


7.5x10" 


-17 


1, 10 


13.0 


4.2x10" 


-11 


0.14 


7.3x10" 


-17 
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